Quantum Dissipative Dynamics of Entanglement in the Spin-Boson Model 



O 
O 

(N 

-«-> 
o 

O 



Jong Soo Lim, 1 Mahn-Soo Choi, 2 ' 3 -[] and M. Y. Choi 1 - 4 

1 Department of Physics, Seoul National University, Seoul 151-747, Korea 
Department of Physics, Korea University, Seoul 136-701, Korea 
3 Department of Physics and Astronomy, University of Basel, 4056 Basel, Switzerland 
J 'Korea Institute for Advanced Study, Seoul 130-722, Korea 

We study quantum dissipative dynamics of entanglement in the spin-boson model, described 
by the generalized master equation. We consider the two opposite limits of pure-dephasing and 
relaxation models, measuring the degree of entanglement with the concurrence. When the Marko- 
vian approximation is employed, entanglement is shown to decay exponentially in both cases. On 
the other hand, non-Markovian contributions alter the analytic structure of the master equation, 
resulting in logarithmic decay in the pure dephasing model. 
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The dynamical properties of an open quantum system 
are conveniently characterized by its decoherence times, 
viz, how long the system can maintain coherence in a 
given superposition state. It is customary and useful to 
distinguish further two limiting decoherence processes, 
relaxation and dephasing, although realistic decoherence 
phenomena appear with both processes admixed. In 
many cases of quantum information processing, a quan- 
tum device consists of two-level systems, for which rates 
of the two processes are designated by 1/Ti and I/T2, 
respectively, in the convention of nuclear magnetic res- 
onance. When the system can be partitioned into two 
parts, one can exploit another interesting property of 
the system, the entanglement between the two parts. 
As highlighted by the EPR paradox [2| and Bell's theo- 
rem 0, |j| , entanglement has been regarded as one of the 
most fundamental properties of quantum physics. In re- 
cent decades, entanglement has attracted renewed inter- 
est as a key resource for quantum information processing 
such as quantum dense coding p|, quantum teleporta- 
tion [jj, and quantum computation Since the first 
experimental generation of the entangled photon pairs 
0, a great number of researches have been devoted to 
the reliable creation and maintenance of the entangle- 
ment, e.g., between far-separated photons 0, between 
solid-state qubits [stllfj). and beween macroscopic atomic 
ensembles (llj . 

A natural question is then how the entanglement be- 
tween the two partitions evolves in time. In particular, 
it will be useful to examine if such dynamics of the en- 
tanglement can be characterized by a separate time scale 
and if so, how it is related to the decoherence times (e.g., 
Xi and T 2 ) of the system. In this paper, we investi- 
gate the dissipative dynamics of two "spins" (two-level 
systems) coupled to environment within the spin-Boson 
model both in the Markovian and non-Markovian limits, 
and examine the characteristic time scale over which the 
entanglement decays. It is found that entanglement may 



decay exponentially or logarithmically, depending on the 
presence of spin relaxation. 

We consider two spins, coupled respecti vely to sepa- 
rate ("local") baths of harmonic oscillators [l^- Follow- 
ing Ref. Il'il we adopt the spin-boson Hamiltonian 



H = H S + H l 
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SB, 



with 
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describing isolated spins, 
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baths of harmonic oscillators, and 
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coupling between the spins and the oscillators, where 
(j = 1,2 and fi = x,y, z) are Pauli matrices representing 
the jth spin and b^ m /bj m are creation/annihilation op- 
erators for the oscillator of mode m coupled to the jth 
spin. The unit system fi — fc^ = 1 is used through- 
out this work. The dynamics of the spins at time scales 
longer than the bath correlation time lo~ x , which is of 
our concern, does not depend on the details of the cou- 
pling constants gj, m an d the oscillator frequencies u)j, m - 
It suffices to characterize each bath collectively in terms 
of the bath spectral density function |l3j 



J jM = E \93,m\ 2 S(uj - UJ hm ). 
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In this work, we focus on the Ohmic model of dissipation 
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where ctj is the dimensionless damping parameter of the 
bath coupled to the jth spin. 

In spite of the simple form of the Hamiltonian, the 
dynamics of the spin-Boson model is highly non-trivial 
and has been the subject of a number of works [l3[ . 
Note in particular that the two spins in the system are 
completely decoupled with each other, each interacting 
separately with its local bath. In equilibrium, the den- 
sity matrix ps of the two spins is therefore separable: 
Ps{t — » oo) = pi(oo) (g>p2(oo)- However, starting from an 
initial state ps(0) with a finite amount of entanglement, 
the system displays non-trivial time evolution of the en- 
tanglement contents in ps(t). [See Refs. 0, H, IU E3 for 
initial preparation of entanglement between two spins.] 

To proceed, we diagonalize Eq. J2J), and obtain 



Hsb 



\ E E rt ( 7 ) 

3=1,2 

(?j,mb] >m + gl m bj,m) , (8) 
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+r? sin 



where E 



= tan 



and 



sis. 



j ■ \ •, A;. «, tan-^AjVej), ., 

= x, y, z) are the Pauli matrices in the rotated ba- 
Here we probe the evolution of entanglement in the 
two limiting cases: When [e^ | 3> |Ay|, i.e, Oj rj 0, each 
spin is not allowed to flip its direction: [tj,H} — 0. In 
the analogy to a particle in a double- well potential [l3| . 
this corresponds to the case of an infinitely high potential 
barrier. The presence of the bath thus cannot change the 
population of spin states, merely breaking the coherent 
superposition of different spin states. In this sense we 
call this limit the "pure dephasing" model |14|. In the 
opposite limit \ej\ <C |Aj|, i.e, 6j ~ 7r/2, with the rotat- 
ing wave approximation employed, coupling to the bath 
leads to relaxation of spin polarizations; this corresponds 
to the strong bias potential in the analogy to the particle 
in a double-well potential. We thus call this limit the 
"relaxation" model. 

The evolution of the density matrix p of the total (iso- 
lated) system is unitary and follows the von Neumann 
equation 



dp 
~dt 



— = -i[H,p] = Cp, 



(9) 



with the Liouville operator C — Cs + £b + £-sb corre- 
sponding to the three parts of H. The reduced density 
matrix ps — Tr^p of spins, given by the trace of p over 
the bath variables J3, is not unitary due to the coupling 
to baths. Its equation of motion, obtained from Eq. © 
by taking the trace, reads 



dps 
dt 



£sPs + Vp s 



(10) 



with the dissipative part formally expressed as [if 



VPS® = - dt'TiB \h sb , e QC ^ [H SB , Ps {t') ® p B ) 
Jo 1 

(11) 

where we have introduced the projection operators V p = 
{Ttbp)®Pb and Q = l—V with the equilibrium bath den- 
sity matrix ps = e~ Hs / T /Ttbp^ Hb ^ t at temperature T 
|l5| . We have also assumed p(0) = Ps(0)®Pb for the ini- 
tial configuration . Further, to get analytic expressions 
for ps (t) , we take the first-order Born approximation and 
make the replacement in Eq. I jllj l 



,QC(t-t' 



,Q{C s +C B ){t-t') 



(12) 



which is valid for sufficiently weak spin-bath coupling. 

We measure the entanglement degree in the mixed 
state ps by means of the concurrence 



C{ps) = max \/Ai 



A 4 , , (13) 



where Ai > A2 > A3 > A4 are the eigenvalues of the ma- 
trix Ps(t v ®Ty)p* s {T y ®T y ). The concurrence vanishes for 
a separable state and becomes unity for a maximally en- 
tangled state. While it is possible to investigate C(ps(t)) 
for a general initial state ps(0), we shall usually assume, 
for clearer physical interpretation, a particular class of 
initial states, namely the Werner states: 



p s (0) = W(r) = r|#+)(# 



1 - r 



I, 



(14) 



where |$+) = ® \l) + ||) |T)) is one of the Bell 

states and / denotes the 4x4 identity matrix. The pa- 
rameter r in the range [0, 1] interpolates between the fully 
random state I/A (for r = 0) and the maximally entan- 
gled state |$+) ($+| (for r = 1). 

The master equation for ps, given by Eqs. 1|10[) - (|12[1 . 
takes the form of an integro-differential equation, and 
the state ps(t) at time t depends on the history ps(t') 
with t' < t. We first discuss the Markovian limit, where 
bath correlations decay sufficiently fast compared with 
the relaxation time of the spins. Mathematically, it cor- 
responds to replacing ps(t') in the integrand of Eq. I|ll|) 
by ps(t). With this Markovian approximation, the mas- 
ter equation l|10|) reduces to a linear differential equa- 
tion with time-independent coefficients. Then expected 
is exponential behavior of ps(t) and accordingly of the 
concurrence C[ps(t)], from which one can extract the 
characteristic time scale of the entanglement dynamics. 

For the pure dephasing model (Aj = 0), the dissipative 
part given by Eq. i|ll|) reads 



ZV>s(i) = y> [T}ps(t)T*-p s (tj\ 



(15) 



where jj = 2-na.jT is the dephasing rate of the jth spin. 
Given an initial state ps(0) = W(r), it is straightforward 
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to find the solution ps(t) of Eqs. (|10|) and l|15|) . It leads 
to the concurrence 



Cu.l) <->(r-~ 



-4(71+72)* 



1 - r 



(16) 



where 0(x) is the unit step function. In particular, for 
the maximally entangled initial state (r = 1), the con- 
currence becomes 



which decays exponentially with the rate proportional to 
the dephasing rate 71 + 72- Notice that the concurrence 
does not decay at zero temperature, where the dephasing 
rates vanishes. 



C(r=l,t) 



-4(71+72)* 



(17) 



For the relaxation model, Eq. (flip reduces to 
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where 7j = 2tt J j(Ej) is the dephasing rate, 

ME) 
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8Ei = Pr 



dE coth 



E 
2T 



E 4 - E 



(19) 



is the environment-induced level shift [15| . and ro(-E) = 
[ e s/T _ I]- 1 i s the Bose distribution function. It is again 
straightforward to find the solution ps(t) of Eqs. H10|) 
and Ijl8|) . which gives the concurrence (for 71 = 72 = 7 
and T = 0) 



C(r,t) = 0(r- l/3)e 



-27* 



e 7 *r 



=7* 



1 + r 



(20) 



It is observed that the concurrence again decays exponen- 
tially and the rate is determined by the dephasing time 
of a separate spin. For the maximally entangled initial 
state, it is simply given by 

-27* 



C(r=l,t) = e 



(21) 



We now go beyond the Markovian approximation (still 
within the first Born approximation). In the non- 
Markovian case, it is convenient to take the Laplace 
transform 



Ps{z) 



dt e- zt p s (t) 



(22) 



and solve the resulting algebraic equations in the place 
of Eqs. (|10fl and (|llfl . Then the dissipative dynamics of 
the system is determined by the pole structure of ps(z) 
on the complex z-plane 16]: Whereas an isolated simple 
pole at z = Zk with residue rk/2?ti contributes an expo- 
nential part rfee Zfctfc to the evolution of ps(z), a branch 
cut of p(z) produces non-exponential behavior. 

The dynamics of ps{t) in the non-Markovian limit is 
governed by the correlation function 



c i (t)= / dwJjM 

Jo 



colli ( -^pj cos(wi) + i sin(cji) 



(23) 
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FIG. 1: Contour for the dephasing model. The branch cut 
has been chosen for computational convenience. 



of the quantum Langevin force ^[j^^ + g* j m b^ m ]. 
In particular, the solution ps{z) involves the respective 
Laplace transforms Cj (z) and C"{z) of the real and imag- 
inary parts of Cj (t) , which are given explicitly by 



cos(z/lo c ) + sin(z/w c ) C\{z/lo c ) 
sin(z/w c ) — cos(z/w c ) Si(z/u c ) 

(24) 

Here Si and Ci are the sine- and the cosine-integral, re- 
spectively. 

For the pure dephasing model, the solution p(z) in- 
volves only C'j(z). Note that the cosine integral Ci(z) in 
Eq. (|24() behaves as In z in the limit z — > while the sine 
integral is analytic. This leads to a branch cut with no 
isolated simple poles on the complex z plane, as shown 
in Fig. ^ K is remarkable that the solution has extra 
structure not present in the Markovian approximation. 
In particular, the absence of isolated poles at T = indi- 
cates that there is no exponential decay in ps(t) and the 
branch cut contribution gives the leading correction. 
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Figure [21 shows the time evolution of the concurrence 
for several different initial states (r = 1/2, 3/4, and 1). 
The intermediate and long time behaviors are given by 



lnw c t and (\nui c t) , respectively. Specifically, for the 
maximally entangled initial state, we have 



C (r,t) = — 



1 

7TO 



(4a + l)7r - 2 tan" 



1 - 2a (7 - m(cj c t)) 



2na 



(4a - 1)tt + 2tan~ 



1 -2a (7- ln(u; c t))' 



27ra 




ln(w c t) 



FIG. 2: Time evolution of the concurrence C(r, i) for the pure 
dephasing model. 
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FIG. 3: Contour for the relaxation model. 



Note that unlike the Markovian case, the concurrence 
decays, albeit slowly, even at zero temperature. 



For the relaxation model, the solution ps{z) involves 
Cj(z) as well as Cj(z). This leads to two branch cuts on 
the complex z plane, as shown in Fig. [3J Another differ- 
ence from the pure dephasing model is the existence of an 
isolated pole. A closer investigation reveals that in the 
long-time limit contributions from the branch cuts are 
negligibly small compared with those from the isolated 
pole [16J. Accordingly, ps(t) in the time domain decays 
exponentially, and the behavior of the concurrence is sim- 
ilar to that in the Markovian limit. 



In conclusion, we have used the generalized master 
equation to study quantum dissipative dynamics of en- 
tanglement in the spin-boson model. The pure dephas- 
ing model and the relaxation one have been considered, 
and entanglement for both models has been shown to 
decay exponentially in the Markovian limit. When non- 
Markovian contributions are taken into account, on the 
other hand, logarithmic decay has been revealed for the 
dephasing model. 
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